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ABSTRACT. In this short note we comment about some criticisms - appeared in 
t^- ' a recent paper by Iguchi et al - to our previous works on gravitational collapse of 

perfect fluids. We show that those criticisms are incorrect on their own. 

o 

Gravitational collapse of perfect fluid is a widely treated topic, especially with 
numerical simulations. To perform an analytical study, studying the line element 
in comoving coordinates turns out to be a hard task, and the so-called area-radius 
coordinates reveal extremely useful for this purpose, though on the other side less 
intuitive for understanding the underlying physical situation. In the following we 
are going to comment about some criticisms to our work EEI appeared in a recent 
. post 1 3 1 by Iguchi et al. 

Dealing with a perfect fluid with linear equation of state p = ae (let us take 
a G (0, 1)), yields the following expression for the matter density p: 
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where \I/(r, R) is the Misner-Sharpe mass function. In our works we basically rely 
on an analytical form for the mass function ^(r, R), that we assume to be Taylor 



^jy expandable at (0, 0): 

*(r,R) = |(r 3 - -^V) + £ *«t'# + ... 
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where the coefficients of third order are found out from the relation, involving R' - 
the derivative of R with respect to the comoving variable r - 
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using the initial position R = r at the first comoving instant t = 0. 

Iguchi et al claim that, with such an expansion for the mass term, the central 
density is constant. But of course one has 
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and the sum above may contain, for instance, terms like r 1 / R 2 , which are obviously 
not regular passing to the limit as (r, R) goes to (0,0). The above expression, 
moreover, shows that the energy of non central shells (i.e. fixing r = r > 0) does 
diverge, if terms with negative power of R axe present. 

On the other side, one may argue that, fixing r = above, terms with negative 
power of R disappear. But this line of reasoning doesn't show to be correct, since 
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r = axes of the (r, R) plane is not a physically meaningful direction in computing 
the limit of any quantity as (r, R) goes to (0, 0). Indeed, passing from comoving 
system (t, r) to (r, R) plane, the condition R(t, 0) = maps the set {(£, 0)} into a 
single point - the origin of (r, R) system, namely. Fixed a shell r , the behavior of 
some quantity (energy density, say) can be studied for any R e [0, r ]. Therefore, 
the only point in the central shell r = which makes physically sense is it! = 0, 
where obviously all quantities considered in (r, R) coordinates have no meaning. 
In short, all the relations involving coordinate changes are valid only where this 
coordinate change in regular - then r > in our case - and area-radius system 
origin (r = 0, R = 0) cannot be considered either a regular or a singular point. 

Again, Iguchi et al claim that the equation for R' above implies, with an analyt- 
ical mass function, that R' goes like (r / R) 2 for all comoving time t. One can see, 
instead, that the limit to (0, 0) depends on the function ^ :R /R 2 , whose expression 
we derived above, and which is not regular near the origin. By the way, one can 
easily check that also marginally bound dust solutions in area-radius coordinates 
leads to a similar paradox: indeed, in this case it is 

^ r T / r \ 1/2] / r \ V2 

R' = — R-r(-) + f - ) 
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and since |^ = £ one finds R' = (f), from which one should in principle argue 
that R = r. This apparently paradoxical results (which is actually true only up to 
central singularity formation, i.e. Vt < t s (0), as we well know from the analysis in 
comoving coordinates) again comes out since (r, i?)-system is not able to make a 
distinction between the singular and the regular central shell. 

Actually area-radius coordinates, though may be of great help in carrying out 
some of the calculations, are not helpful to determine whether the central shell gets 
singular or not. For this aim, comoving coordinates must be involved again to prove 
that the limit of the singular curve t s (r) is finite (and nonzero of course) as r — > 0. 
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